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The auto-correlation between two members of a galaxy population is symmetric under the inter¬ 
change of the two galaxies being correlated. The cross-correlation between two different types of 
galaxies, separated by a vector r, is not necessarily the same as that for a pair separated by — r. 
Local anisotropies in the two-point cross-correlation function may thus indicate a specific direction 
which when mapped as a function of position trace out a vector field. This vector field can then be 
decomposed into longitudinal and transverse components, and those transverse components written 
as positive- and negative-helicity components. A locally asymmetric cross-correlation of the longi¬ 
tudinal type arises naturally in halo clustering, even with Gaussian initial conditions, and could be 
enhanced with local-type non-Gaussianity. Early-Universe scenarios that introduce a vector field 
may also give rise to such effects. These antisymmetric cross-correlations also provide a new pos¬ 
sibility to seek a preferred cosmic direction correlated with the hemispherical power asymmetry in 
the cosmic microwave background and to seek a preferred location associated with the GMB cold 
spot. New ways to seek cosmic parity breaking are also possible. 

PACS numbers: 98.80.-k 


I. INTRODUCTION 

Considerable evidence has accrued for the past decade 
that large-scale structure in the Universe grew via grav¬ 
itational infall from a nearly scale-invariant spectrum of 
very nearly Gaussian primordial adiabatic density per¬ 
turbations mm- The nature of these perturbations is 
generally accounted for in terms of single-field slow-roll 
(SFSR) inflation [3]. However, even if SFSR is the correct 
explanation, there might be new physics beyond SFSR 
inflation, and there remains the ever-present possibility 
that primordial perturbations are due to something com¬ 
pletely different. 

In order to develop a clearer understanding of the new 
physics responsible for primordial perturbations, we must 
be vigilant in seeking new fossils from the early Universe 
in the form of subtle correlations in the matter distri¬ 
bution. This is the motivation for much of the work on 
non-Gaussianity [3]. This Letter will explore new observ¬ 
ables that can be sought with galaxy surveys. Such work 
is timely given the advent in the forthcoming years of a 
new generation of surveys [SHS] that will map the distri¬ 
bution of galaxies over vast volumes in the Universe. 

Recent work [5] presented a parametrization of the 
most general autocorrelation function—which may de¬ 
pend on the orientation of the two points being corre¬ 
lated as well as their position in space—and showed that 
it could be decomposed into scalar, vector, and tensor 
components. Here we generalize that work to the most 
general two-point cross-correlation function between two 
different galaxy populations. If primordial perturbations 
are Gaussian, then the two-point correlation function 
is statistically isotropic and homogeneous. More gener¬ 
ally, though, the two-point correlation function may, at 
least in some small region of space, be anisotropic. That 


anisotropy may also vary from one point in the Universe 
to another. Ref. [3] parametrized the most general local 
departure from isotropy for an auto-correlation function. 
Since an auto-correlation function for two galaxies sep¬ 
arated by a vector r must be invariant under the inver¬ 
sion r —>■ —r, the most general departure from statistical 
isotropy is parametrized in terms of a symmetric tensor, 
or more generally, in terms of the six degrees of freedom 
that parametrize that tensor. 

If, however, we consider the two-point cross-correlation 
function between two distinct populations, then it is pos¬ 
sible that the correlation function for a galaxy pair sep¬ 
arated by r may differ from that with separation —r. 
There is thus a possibility that the two-point correla¬ 
tion function might “point” in a given direction, and this 
pointing is described by a vector. Our purpose here is to 
point out that evidence of the imprint of a vector field can 
be sought with cross-correlatoins in the galaxy distribu¬ 
tion. Our work differs from that of Refs. unHH], who also 
considered asymmetric galaxy cross-correlations; while 
they considered only those that arise from projection ef¬ 
fects, we consider bona fide asymmetries in the three- 
dimensional distribution. Since observations are done in 
redshift space, though, care must be taken to disentangle 
the 3D asymmetries we consider below from the redshift- 
space asymmetries discussed in Refs. [inHn], and from 
the angular dependence induced by the 3D RSD opera¬ 
tor [l3] . 

Below we describe the parametrization of the most gen¬ 
eral two-point cross-correlation, beginning by reprising 
the parameterization for the auto-correlation function. 
We then discuss the decomposition of these types of local 
departures from statistical isotropy in terms of two trans¬ 
verse and one longitudinal mode. Next, we write down 
the optimal quadratic estimators to be constructed from 
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galaxy surveys to detect these departures from statistical 
isotropy. After that, we discuss several possible physi¬ 
cal mechanisms to generate these types of local depar¬ 
tures from statistical isotropy. We begin this discussion 
with an effect that arises in biased halo clustering, even 
with Gaussian initial conditions, and another that arises 
in halo clustering if there is local-type primordial non- 
Gaussianity. We then describe several early-Universe sce¬ 
narios for such correlations. We point out that a search 
for these anisotropic cross-correlations may be used to 
seek a preferred cosmic direction correlated with the di¬ 
rection indicated by the hemispherical power asymmetry 
[HIIS] or to seek a preferred location associated with the 
GMB cold spot m- Finally, we mention the possibliity 
to construct new probes of parity breaking. 


II. PARAMETRIZATION OF THE 
CROSS-CORRELATION 

We begin by reprising the parametrization of Ref. [5] 
for the two-point auto-correlation function. Suppose we 
have a density held (5(x) from which we construct the 
Fourier components (5(k). The most general two-point 
correlation function for those Fourier components can be 
parametrized as, 

(5(ki)5(k2)) = P(fci)4i,-k. + 

K p 

xh;(K)eA(K)A:lfc^2C3^ (1) 

where <5ki,-k2 is a Kronecker delta, and ^ Kro- 

necker delta that sets ki -|- k 2 -f K = 0. Here p sums 
over the six possible basis tensors for a symmetric tensor, 
and e^j(K) are the six polarization tensors which can be 
written as a trace, a longitudinal mode, two transverse- 
vector modes, and two transverse traceless tensor modes. 
The hp(K) are Fourier amplitudes of the various types 
of perturbations, and fp(ki,k 2 , p) parametrizes the de¬ 
pendence on ki and k 2 (where p is the cosine of the 
angle between ki and k 2 ). The hrst term in Eq. Q in¬ 
dicates the statistically isotropic correlation with power 
spectrum P{k). It was argued in Ref. [3], simply from 
symmetry considerations (and in particular, the depen¬ 
dence of the correlation on the azimuthal angle about the 
direction of K), that terms with the two scalar, the two 
vector, or the two tensor polarizations can arise in infla¬ 
tion if the inflaton is coupled to a new scalar, vector, or 
tensor field, respectively. 

Our purpose here, though, is to consider the additional 
possibility that arises when we cross-correlate two differ¬ 
ent tracers of the matter distribution. The fractional 
density perturbations of these two populations will have 
Fourier amplitudes i5i(ki) and (52(k2), respectively. The 
most general two-point correlations of these Fourier co¬ 
efficients will satisfy a relation like Eq. Q, but the sum 
on p will now be extended from six to nine to account 
for the three possibilities where efj(K) is antisymmetric. 


The three new degrees of freedom in the cross-correlation, 
not allowed with auto-correlations, can be written most 
generally, for ki ^ —k 2 , as, 

(<5i(ki)52(k2)) = E /p (ki, k 2 , h;(K)ep-(ki - k 2 ) ^ 3 . 

K,p 

( 2 ) 

Here, the sum on p runs over the three polarizations 
p = L,x,y, where eL(K) = K and ex,y(K) are two other 
unit vectors orthogonal to K and to each other. Just 
from symmetry considerations, we expect that correla¬ 
tions with p = L could arise if there is some longitudinal- 
vector held—or equivalently, the scalar held from which 
it is derived—coupled to whatever physics determines the 
galaxy distribution. Similarly, distortions with p = x,y 
would require that the galaxy distribution was deter¬ 
mined, at least in part, by some transverse-vector held. 

Eq. ([^ describes, in Fourier space, cross-correlations 
between two different populations, that are antisymmet¬ 
ric in the exchange of the two populations. In conhgura- 
tion space, these correlations trace out a vector held, as 
follows: In any small volume of the Universe, the cross¬ 
correlation could “point” in some given direction; i.e., 
there could be a mean offset between galaxies of type 
I (e.g., more massive galaxies) and galaxies of type 2 
(e.g., less massive galaxies). This thus implies a pre¬ 
ferred direction in that volume. The generality of the 
parametrization in Eq. ([^ allows for the possibility that 
this preferred direction is spatially dependent in such a 
way that global statistical isotropy is preserved on suf- 
hciently large scales. It also allows for the possibility, 
though, for a preferred direction across the entire ob¬ 
servable Universe. This latter possibility clearly requires 
some exotic new physics, but the imprint of a local pre¬ 
ferred direction is not at all exotic, as we will see. 


III. ESTIMATORS 


We now describe how to measure these vectorial cross¬ 
correlations from a galaxy survey. We will consider 
the longitudinal and transverse components separately 
and begin with longitudinal distortions, supposing that 
/l(A:i, ^ 2 ,/i) is specihed. Suppose further that the phys¬ 
ical model specifies that the Fourier components ^^(K) 
are selected from a Gaussian distribution with power 
spectrum Pl{K) = ^lP[{K) of amplitude Al and fidu¬ 
cial K dependence P[{k). Following the same steps as 
in Ref. [5], the optimal estimator for the amplitude 
is given by 


Al=ctIY, 


K 


2[Pl{K)Y 


V 


-1 


hL(K) -PE(K) 


(3) 

and this amplitude is determined with inverse variance, 


-2 


{l/2)Y,[p[{K)f [P£{K)] 

K 


(4) 
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Here, 


MK) 


V fLik,\K-k\,fi)K-{K + 2k) 
Ly ) 2 VPl°\k)P^°WK - k|) 

x^i(k)<52(K-k), (5) 


is the estimator for the amplitude /ii(K), where here 

= k ■ (K — k), and the noise power spectrum for hi(K) 
is 


P£{K) = 


E 


\fL {k,\K - kl fi) k -{KP 2 k) 


1 -1 


2VPl^\k)P^^^{\K-'k\) 


( 6 ) 

In these expressions, V is the volume of the survey, the 
sums are over the wavenumbers k for which the density 
perturbations are measured; Pi{K) and P2{K) are the 
power spectra for the two populations; and is the 
total power spectrum, the sum of the signal and the shot 
noise. If instead we consider distortions of the transverse- 
vector type, then the optimal estimator for the amplitude 
At for the power spectrum for the transverse-vector field 
will be as above, but with K replaced by ei(K) and an 
additional sum over the two transverse polarizations. 


IV. SOURCES OF ASYMMETRIC GALAXY 
CROSS-CORRELATIONS 

A. Biased halo clustering 

We now show that asymmetric cross-correlations arise 
in the standard model for biased halo clustering from 
Gaussian initial conditions. Let i 5 (x) be the primordial 
fractional mass-density perturbation, and suppose that 
this density field is populated by two different types of 
halos. These could be, for example, high-mass halos (e.g., 
those that house giant ellipticals, clusters, etc.) and 
low-mass halos (e.g., those that house Milky-Way type 
galaxies, dwarf galaxies, etc.). Let ni(x) be the frac¬ 
tional perturbation in the number density, as a function 
of position x, of halos of type 1 and similarly for n2(x). 
Since the abundance of halos in some region is a nonlin¬ 
ear function of the local mass density, there are nonlin¬ 
ear relations [T^, 'Ui(x) = &i( 5 (x) -|- ci[( 5 (x)]^ -!-••• and 
n2(x) = 62 < 5 (x) -I- C2[(5(x)]^ -!-••• between the halo and 
mass densities, where bi and &2 are the usual linear-bias 
parameters, and Ci and C2 are nonlinear-bias parameters. 
The nonlinear relation between rii and S{x) implies that 
there will be a three-point correlation, 

(ni(ki)n2(k2)5(K)) = 2P{K) [&2CiP(/c2) + &iC2R(fci)] 

( 7 ) 

If we antisymmetrize in ki and k2, we find an antisym¬ 
metric part, oc 2P{K){b2Ci — biC2) [P{ki) — P{k2)] to 
this three-point function. We now consider the squeezed 


limit, k = ki ^ —k2 and k ~ ki,k2 ^ K, appropriate 
if we are considering the small-scale clustering of halos 
in the presence of a low-pass-filtered density field. The 
part of the bispectrum antisymmetric in ki and k2 is 
then 2 {biC 2 — b2Ci)P{K)[dP{k)/dk]K ■ k. We thus see 
that in the presence of a long-wavelength mode of i 5 (x), 
the high-pass-filtered ni(x) and n2(x) have a Fourier- 
space cross-correlation precisely of the form in Eq. 
with /l = [dP{k)/dk]{biC2 — h2Ci){K/k). The longitudi¬ 
nal, rather than transverse, mode is as expected in the 
presence of the long-wavelength density perturbation. 

Existing measurements [TM^ of nonlinear-bias pa¬ 
rameters from galaxy bispectra have relied en¬ 

tirely on measurement of auto-correlations. Cross¬ 
correlations have also been used to infer bias parameters 
[? 7 ll 30 | . but those works considered only the linear bias. 
If a galaxy sample is broken up into two different popu¬ 
lations, the linear- and nonlinear-bias parameters for the 
two populations can be measured with the bispectra for 
the two different populations. The symmetric part of 
the cross-correlation between these two populations then 
provides a measurement of the combination biC2 P 62C1, 
and the antisymmetric part of this cross-correlation pro¬ 
vides 62C1 — biC2. The latter provides some incremental 
improvement in statistical power, and it can also com¬ 
plement other measurements of the bias parameters and 
be used to check the validity of the biasing-model as¬ 
sumptions behind the analyses. Indeed, from Eq. Q 
and what follows, we see that for a power-law P{k), the 
ratio between the symmetric and asymmetric signals is 
cx (62C1 — biC 2 )K/ 2 {biC 2 P &2Ci)/c, which is generically 
small in the squeezed limit k ^ K. However, for popula¬ 
tions with negative non-linear bias on certain scales (see 
e. g. [33]), it is certainly possible that the bias parame¬ 
ters contrive to yield 61C2 P 62C1 « 0 , in which case the 
asymmetric information becomes particularly beneficial. 


B. Primordial non-Gaussianity 


Asymmetry may arise even with linear biasing if 
there is primordial non-Gaussianity of the local type. 
In this case, the linear-bias parameters for the two 
halo populations may be scale dependent, ni(ki) = 
bi{ki)S{ki) and n2(k2) = b2{k2)S{k2). Here bi{k) « 
bi P 2 {bi — l) 5 cfn\^m'^{a.H/k)'^ /2 with the standard 
scale-independent bias bi and some collapse thresh¬ 
old 5 c [3n 1311 131] • On the other hand, the matter- 
density perturbation is expected to develop a nonzero 
bispectrum B{ki,k2, K) due to nonlinear gravitational 
growth | 55 H 57 ] (this should not be confused with 
the primordial bispectrum). The three-point func¬ 
tion (ni(ki)n2(k2)(5(K)) oc bi{ki)b2{k2)B{ki, k2, K) 
is thus noninvariant if we exchange ki and k2, 
and we again recover an antisymmetric correlation oc 
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[bi{ki)b2{k2) - bi{k2)b2{ki)] B{ki, k2, K)/P{K), with 


h 


3(6i 


^2)^c/nl^m 


{aH)^K B{ki,k2,K) 
“fcffcf P{K) 


( 8 ) 


C. Primordial longitudinal vector field 

We now describe an exotic scenario that could give 
rise to an anisotropic two-point correlation. Consider a 
curvaton-like inflation model in which primordial density 
perturbations are due not to fluctuations in the inflaton, 
but to fluctuations in some spectator field. Suppose fur¬ 
ther that this curvaton-like field is complex and coupled 
to a new 17 ( 1 ) gauge field. There will then be charge- 
density fluctuations that arise during inflation, and if the 
17 ( 1 ) symmetry is broken after inflation, those charge- 
density fluctuations may survive. It is then conceivable 
that dark matter may be like baryonic matter in our Uni¬ 
verse, in which there are light electron-like dark-matter 
particles and heavier proton-like dark-matter particles of 
opposite charge. If so, then galaxies that form in re¬ 
gions with a positive charge-density fluctuation may be 
different than those that form in regions with negative 
charge-density fluctuations. The cross-correlations be¬ 
tween these two different types of galaxies—which for ex¬ 
ample may appear to us as brighter or fainter galaxies— 
should then trace out the large-scale electric-type field 
associated with the early-Universe 17 ( 1 ) gauge field (see 
Fig. [^. Such a longitudinal-type cross-correlation could 
in principle be distinguished from that due to halo clus¬ 
tering from the different dependences on k. Unlike the 
halo-biasing case, in which small-scale directionality in¬ 
dicated by anisotropic cross-correlations are correlated 
with the long-wavelength density field, these anisotropic 
correlations would have no cross-correlation with the 
large-scale density field. 


D. Primordial transverse vector field 


K of the vector field to two scalar-field modes of wavevec- 
tors ki and k2 would be symmetric under ki k2, one 
could again construct a model with two dark-matter com¬ 
ponents, sourcing two different galaxy populations in the 
respective dark matter halos made up from each compo¬ 
nent, in a way that makes the cross-correlation in one 
of the transverse directions asymmetric (see illustration 
in Fig.[^. These transverse-vector cross-correlations will 
be straightforward to measure with a galaxy survey and 
any non-null result would be extraordinarily interesting. 
The measurement will also be useful as a unique null test 
for systematic effects. 


E. The CMB power asymmetry and cold spot 

There is evidence for a hemispherical power asymmetry 
in the cosmic microwave background [Ml [T^ and it is 
important to investigate whether this preferred direction 
shows up anywhere else. If the CMB power asymmetry 
is due to the long-wavelength modulation of primordial 
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What about anisotropic cross-correlations of the trans¬ 
verse variety? As we have seen above, modulation 
of small-scale galaxy-density fluctuations by a long- 
wavelength density mode arise naturally, but those are 
of the longitudinal variety. Any small-scale anisotropic 
cross-correlation that traced out a transverse-vector field 
would thus provide clear indication of new physics. There 
are indeed a number of inflationary scenarios that involve 
new vector fields |38j and models of gravity in which the 
vector degrees of freedom in the metric are brought to life 
[89j , and it is reasonable to surmise that one of these pro¬ 
vide an imprint of the type we consider. To see how this 
may work, suppose the inflaton (or curvaton) (j) is coupled 
to a vector field through a term A‘'\di_i(j)){d^(j)), 
where the parentheses in the superscripts indicate sym- 
metrization (antisymmetrization would make this term 
vanish). Although the coupling of a mode of wavenumber 


0© 0© 0© 

-► 


FIG. 1: We consider models with two components of dark 
matter with some dipole interaction. Top: An “electric field” 
in the longitudinal direction leads to the arrangement in the 
pattern shown above. Then at the separation scale corre¬ 
sponding to the distance between and “-I-”, there is an 
asymmetric cross-correlation between the and “-I-” pop¬ 
ulations, along this longitudinal direction. The estimators 
derived in this work are designed to detect such a transverse 
vector asymmetry. Botom: A different setup, for which the 
asymmetry has a sinusoidal modulation in the transverse di¬ 
rection, leading to a transverse asymmetric cross-correlation. 
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perturbations |40H45j . then anisotropic cross-correlation 
of the type discussed above for biased halo clustering 
should exist, but with just one long-wavelength mode 
/ii(K) of wavenumber K —>■ 0 in the direction of the 
CMB power asymmetry. In simple terms, this Fourier 
amplitude hi(K) seeks whether a less massive galaxy is 
more likely to be found on one side of a more massive 
galaxy than on the other side. 

Measurements from Planck dsj m] also conhrm the 
existence of a cold spot in the CMB m- An asymmetric 
cross-correlation that traces out a hedgehog configuration 
surrounding the cold spot could conceivably arise either 
from a large void [4(i) or from exotic physics m that may 
be responsible for the cold spot. A template for such a 
configuration could be constructed from the estimators 
/il(K) discussed above. 


F. Parity breaking? 


There are also novel probes of cosmic parity breaking 
that can be constructed. Estimators for the two trans¬ 
verse modes /ip(K) for each wavevector K can be con¬ 
structed analogously to Eq. ([^ by replacing the vector 
K that appears in the dot product with one of the two 
transverse polarization vectors. Those two linear polar¬ 
ization vectors ix and iy can be replaced by circular po¬ 
larizations e±(K) = 2“^/^[ea;(K) ± ley(K)]. Estimators 
for the Fourier amplitudes of these helicity states are then 
given by 


AK) = pE{k)J2 


2EP‘°‘(fc)P‘°*(|K-k|) 


X(5i(k),52(K-k). 


It is then straightforward, following the power-spectrum- 
amplitude estimators discussed above, to construct esti¬ 
mators to check for an asymmetry between the power in 


right and left circularly polarized modes. Again, while a 
null result would not be too surprising, a positive result, 
if found, would be revolutionary. 


V. CONCLUSIONS 


We have shown that cross-correlations between differ¬ 
ent galaxy populations may “point” in a given direction. 
The vector field from this pointing can be reconstructed 
from galaxy clustering, and it can be decomposed into 
longitudinal and transverse components, and also into 
components of positive and negative helicity. We dis¬ 
cussed several physical mechanisms that may give rise to 
these departures from local statistical isotropy. These in¬ 
clude nonlinear halo biasing, an effect that should arise 
in the standard model of halo clustering, and local-type 
non-Gaussian primordial perturbations. Other possibli- 
ties include couplings of the primordial field responsible 
for primordial perturbations to a new vector field, as well 
as a coupling to some field that might account for the 
hemispherical power asymmetry or cold spot in the CMB. 
The clustering effects from halo biasing constitute con¬ 
crete predictions of the standard halo-clustering model 
that can be measured and used to test that standard 
model. Although the exotic possibilities we discussed are 
long shots, they will be easily sought in forthcoming sur¬ 
veys and would, if discovered, be quite remarkable. 
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